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ABSTRACT 

Recent results have suggested that the density of baryons in the Universe, Ob, is 
much more uncertain than previously thought, and may be significantly higher. We 
demonstrate that a higher J^b increases the viability of critical-density cold dark mat- 
ter (CDM) models. High baryon fraction offers the twin benefits of boosting the first 
peak in the microwave anisotropy power spectrum and of suppressing short-scale power 
in the matter power spectrum. These enable viable CDM models to have a larger Hub- 
ble constant than otherwise possible. We carry out a general exploration of high JIb 
CDM models, varying the Hubble constant h and the spectral index n. We confront 
a variety of observational constraints and discuss specific predictions. Although some 
observational evidence may favour baryon fractions as high as 20 per cent, we find 
that values around 10 to 15 per cent provide a reasonable fit to a wide range of data. 
We suggest that models with fiB in this range, with h ~ 0.5 and n ~ 0.8, are currently 
the best critical-density CDM models. 

Key words: large-scale structure of the Universe ~ cosmology: theory - dark matter 



1 INTRODUCTION 

There is presently greater uncertainty regarding the baryon 
density of the Universe than there has been for many years. 
The standard theory of big bang nucleosynthesis (BBN), 
developed over more than 40 years and applied exten- 
sively to local abundance data, has long been considered 
as fixing the baryon density to a satisfyingly high level 
of accuracy. A typical quoted value a few years ago was 
n^th? = 0.0125 ± 0.0025, where the uncertainty is t o be in- 
terpreted a s something like 95 per cent confid ence (Walkei 



et al. 1991; Smith, Kawano & Malaney 1993). Here 57b is 



the baryon density in units of the critical density, and h is 
the present Hubble parameter in units of 100 kms~^ Mpc~^. 
However, recently it has been acknowledged that there are 
a range of possible systematic unce rtainties in determina- 
tions of the primordial abu n dances ( Wjlson_fc_jtood_j_994[ 
Sassel ov fc Goldwirth 19951 ; pkillman, Terlevich fc Garnet] 



Added to that, it is now becoming possible to mea- 
sure element abundances at high redshift as well as lo- 
cally, through absorption systems in the spectra of distant 
quasars. Results for deuterium have so far been inconclusive 
(Songaila et al. 1994; Carswell et al. 1994, 1996; Wampler 
et al. 1996; Rugers & Hogan 1996a,b; Tytler, Fan & Buries 
1996; Buries & Tytler 1996), with some very high abun- 
dances reported (implying a low baryon density) and also 
some low ones. There may be some reason to favour the 
lower determinations, in that most of the likely systematic 
effects, such as interloper hydrogen clouds, bias the abun- 
dance estimates upwards. However, only determinations in a 
large number of separate systems will truly resolve the issue. 
The mo st recent determination giving a low deuterium abun - 



dance (Tytler, Fan fc Buries 1996; Buries fc Tytler 1996) 



1995 



Scully et al. 1996), which have led to a broadening 



of the preferred interval, particularly towards the upper 
end. Recent quoted ranges, at 95 per cent confidence, in- 
clude 0.010 < Q.Bh^ < 0.022 by Copi, Schramm fc T urner 
(1995a,b), with a similar result by Turner et al. (1996), and 
0.0125 < SIb/i^ < 0.0275 by Hata et al. (1996). Kernan 



& Sarkar (1996) quote a 95 per cent confidence uppe r limi t 
Q.B.h? < 0.032, while an analysis by Krauss & Kernan (1991:) 
has also suggested higher baryon densities, although without 
quoting a specific range. 



yields the resuh ^b/i = 0.024 ± 0.002 ± 0.002 ± 0.001, with 
the lo uncertainties being statistical, systematic, and the- 
oretical respectively. This is still consistent with the more 
recent determinations using local abundances, while rein- 
forcing the view that the baryon density may be higher than 
previously thought. 

Within the conventional picture, the new deuterium re- 
sults can only be made consistent with BBN if the primor- 
dial helium abundance Up has been underestimated; Q.B,h? — 
0.025-0.035 corresponds to 25 per cent helium, rather than 
the usual 23-24 per cent. However, this seems to be entirely 
within the bounds of possibility (e.g. Olive & Steigman 1995; 
Sasselov & Goldwirth 1995; Buries & Tytler 1996), given 
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possible systematic effects. The low deuterium values also 
imply very little destruction of primordial deuterium, in or- 
der to be consistent with the interstellar me dium values. 



whi ch give a constraint around Qsh < 0.031 ( McCuUough 



1992 



nsky et al. 1995 ). Hence values of 57b <; 0.15 (assum- 



ing h ~ 0.5) begin to run into quite firm limits. However, it 
may also be possible to relax the nucleosynthesis bound on 
Qb with new particle physics, an example bein g the decay- 
ing tan neutrino proposal of Gyuk & Turner (1994), or by 
allowing inhomogeneities in the baryon-to-photon ratio (see 
e.g. Mathews, Kajino & Orito 1996 and references therein). 

The other salient observational issue is the question 
of th e baryon fraction in clus t ers of galaxies (IWhite et al. 
1993bl IWhite fc Fabian 1995|; |Elbaz, Arnaud fc Bohringer 



1995 



i 



arkevitch et al. 199f: ), which consistently gives values 



which are high for a critical-density unive rse. F or example, 
a recent compilation by White & Fabian (1995) gives 



^ ^ 0.14+" °* 



-0.04 



(4) 



-3/2 



(1) 



again at the 95 per cent level. There does seem to be vari- 
ation in this quantity for individual clusters, but it is still 
unclear whether this is due to systematic effects. Certainly 
some clusters appear to have a baryon fraction as high as 
20 per cent (see e.g. Mushotzky 1995). At the moment there 
is no obviously reliable lower limit to the cluster baryon 
fraction, although the values adopted by Steigman & Felten 
(1995), Qb/^o > 0.2(/i/0.5)"^/^ and by Evrard, Metzler 
& Navarro (1996), Qb/^Io > 0.11 (/l/0.5)"^/^ are typical of 
the range. Assuming there has been no significant segrega- 
tion of the baryons, and that there has been no serious sys- 
tematic underestimation of the mass of galaxy clusters (see 
e.g. Gunn & Thomas 1995; Balland & Blanchard 1996), this 
is hard to reconcile with a critical-density universe, unless 
f^B has a higher value than given by the canonical BBN 
numbers. 

There are some additional considerations which may 
favour an increase in the baryon fraction, although they 
would not make strong arguments on their own. One ex- 
ample is that a higher value of Qb may bring simulations 
of Lyman alpha absorption systems more in line with data 



Katz et al. 1996 



Miralde-Escude et 



( Hernquist et al. 1996 

al. 199^). Another example is the apparently high baryonic 
mass fr action in the halo of our Galaxy, in the form of MA- 



CHOs (A,Icock et al. 1996; Bennett et al. 1996b). And a 



third is that higher 57b may help proto-galaxy disks be more 
self-gravitating, thereby alleviating a problem with over- 
concentration of baryons in simulations of galaxy formation 
(see e.g. Steinmetz 1996). 

We are interested in the implications of these changing 
notions about fie for models of structure formation , and in 



part icular for the cold dark matter (CDM) rn odel ( Peeble; 
1982|; |B|Lumenthal et al. 1984 [Pavis et al. 1985|) in the case of 



a critical-density universe (see Davis et al. 1992a and Dodel- 
son. Gates & Turner 1996 for recent reviews) . We will argue 
that a factor of two increase in the predicted baryon content 
from its circa 1993 value has a far from insignificant effect 
on the viability of such models, and indeed that there exist 
models from this class which provide a reasonable fit to the 
current data. The bottom line will be that models which are 
selected to conform with a wide range of constraints have: 



(1) moderately tilted initial conditions; (b) h ~ 0.5; and (c) 
f^B centred around a value of say 12 per cent. 

The outline of the paper is as follows. In the next sec- 
tion we give an overview of the situation regarding critical- 
density CDM models, and in the following Section define our 
terminology and motivate our choice of model parameters. 
We then compare these high Ob models with the available 
data on structure formation and on the cosmic microwave 
background (CMB) in Sections 4 and 5 respectively. We con- 
clude with a discussion of some of the strengths and weak- 
nesses of these models and future key observations. 



2 THE BIG PICTURE 

In order to obtain a critical-density universe of sufficient 
age, it is necessary to choose a Hubble parameter not much 
greater than h = 0.5 (for example, an age greater than 
12 Gyr requires h < 0.55). Similar values for h are also 
necessary if one is to obtain a satisfactory structure forma- 
tion model ( [White et al. 1995b| ; [Liddle et al. 1996b| ). Some 
recent observational analyses have l ent support to such a 



low Hubble consta nt (Schaefer 1996 



Tammann et al. 199f 



Branch et al. 1996), and such values e ven r emain marginally 
consistent with the Freedman et al. (1994) measurement of 
h — 0.80±0.17. With h around 0.5, the predicted age of the 
uiuverse appears to be consistent with the ages of the oldest 
globular clusters (Bolte & Hogan 1995; Chaboyer et al. 1995, 
1996; Jiminez et al. 1996). For this value of h, the old nu- 
cleosynthesis determinations would indicate a universe with 
about 5 per cent baryons, while the higher values favoured 
more recently would be more like 10 per cent. As we shall 
see, this increase is far from insignificant, and indeed could 
be regarded as a new 'standard' value^. The cluster num- 
ber, centred around 15 per cent and going up to above 20 
per cent, suggests the possibility of an even higher value. 

If one is determined to retain a critical-density universe, 
it is extremely interesting that there is now strong evidence, 
from a structure formation perspective, in favour of these 
higher f^B values. The demonstration of why this makes 
such a difference is the principal goal of the present pa- 
per. The gist of the argument is as follows. It is well known 
that, despite the strong appeal of simplicity and physical 
motivation, the standard CDM model (sCDM), based on a 
scale-invariant initial spectrum (n = 1), a Hubble constant 
h = 0.5 and with 5 per cent baryons, is unable to fit all 
the observational data. This is largely because, when nor- 
malized to reproduce the microwave anisotropics det ected by 
the Cosmic B ackground Explorer (COBE) satellite (Bennett 
et al. 1996a), the power spectrum has too much short-scale 



power. If one is to retain the CDM hypothesis, modifications 
are required, and two have received particular attention, the 



first being to reduce the Hubble parameter (Bartlett et al 



1995) an d the second to tilt the spectrum of density pertur- 
bations (Vittorio, Matarresse fc Lucchin 1988t [Bond 1992 



Liddle, 



Lyth 



cia et al. 1993 



Sutherland 1992[ |Cen et al. 1992| ; [Muciac-, 
Starobinsky 19921 ). The former 



Polarski 



has the drawback that it appears necessary to reduce h to 



In fact, many A'^-body simulations and other CDM calculations 
of the 1980s adopted this round number as a working value, and 
it is still sometimes used in this context. 



© 1996 RAS, MNRAS 000, 000-000 



CDM models with high baryon content 3 



about ( .35 (Liddle et al. 1996b), which is lower than most 1984 ; Kofman fc Starobinsky 1985| [Efstathiou, Sutherland 



researcl lers are willing to accept. The latter has the prob- 
lem the t_^f_OTie_^ilts_^he_s£ectrui^^uffi£^^ 
unwantp^_shprt^scak_^o^^^^^hfi_jn^te£_^ 



: Maddox 199C SuKinohara fc Suto 1991 ; Efstathiou, Bond 
fc White igjllKofman, Gnedin fc Bahcain^"!; iKrauss & 
Turner 1995|; lOstriker fc Steinhardt 1995| ; [Stompor, Gorski 



the tilt is also enough to remove the peak from the CMB fc Banday 1995 



Liddle et al. 



anisotropy spectrum (White et al. 1995b; White 1996). 

A range of CMB data now implies the existence of a 
power-spectrum peak at sub-degree scales (Scott, Silk & 
White 1995; Kogut & Hinshaw 1996), which seems even 
more co mpelling with the recen t Saskatoon-95 (SK95) obser- 
vations f Mettcrficld et al. 1996 ) . This means that purely tilt- 



1996(j ). 



This has been seen as 
the easiest way to make CDM models consistent with the ob- 
servational constraints, as indicated by the large number of 
papers exploring this possibility. We note though that pres- 
sure has been exerted on both the open and flat cosmological 
constant do minated approaches rec ently. The c ombination 



ing the spectrum, while retaining the other standard CDM 
parameters, is beginning to look untenable. The situation 
improves somewhat if one c ombines a smaller t i lt with a 
more modest reduction in h (White et al. 1995b; Liddle et 



of th e SK95 (Netterfield et al. 1996) and CAT (Scott et al 



1996) CMB observations appears to fix the location of the 
acoustic peak at around £ ~ 200, suggesting a flat geometry 
rather than a hyperbolic one. Furthermore, the somewhat 



al. igglblT ^i^^^^^^h^^c^^hi^^^^£^renU^^til^wt 



lower norm a lization of the four-year COBE data ( Bennett 



enough if the currently favoured peak height is confirmed 
by future experiments. 

The most natural way to try to resolve this conflict, 
while staying within the CDM paradigm, is to raise the 
baryon fraction. This has twin benefits. Regarding the ra- 
diation power spectrum, the effect of the extra baryons is 
to enhance the first peak. With enough baryons, this can 
go a long way towards compensating for the loss of height 
introdu ced bv the tilt (the tilt in turn being required if one 



[Gorski et al. 



1996 



Hinshaw 



et al. 1996£|: [Banday et al. 1996[ 

et al. 1996 i" ]Wrlghtet'aL'T996^ xacerbates the problem that 
low-density open universes have a low normalization. Mean- 
while, meas ures of the deceleration parameter using Type la 
supernovae ( Perlmutter et al. 1996| ) suggest a universe which 



is decelerating rather than accelerating, which for a flat uni- 
verse requir es flp > 2/3, a greeing with the quasar lensing 
constraint (Kochanek 1996). In addition, types of observa- 



is to pr pvent the required h being too small), and hence re- 
tain compatibility with the recent CMB observations. The 
second beneflt is that the barvons themselves help suppress 



tions which have historically favoure d high dens i ty, such as 
velocity power spectrum measur es (Dekel 1994; Strauss & 
Willick 199E ; Zaroubi et al. 1996 ) , redshift-space distortions 
(Cole, Fisher & Weinb erg 1995; Taylor fc Hamihon 1996) 



the sho lrt-scale power in the matter power spectrum, since 1994| ), and the absenc e of a prominent break in the non - 



the p roperties of voids ( Nusser fc Dekel 1993 ; Dekel fc Ree; 



unlike the CDM they are unable to collapse until after de- 
coupling. This means that the correct power spectrum shape 
can be obtained for a larger Hubble parameter, or a more 
modest tilt (cf. White et al. 1995b). 

Before embarking on our discussion of critical-density 
CDM models, let us comment on some of the alternatives. 



component of the dark matter is hot ( Shafi fc Stecker 1984J 


Bonometto fc Valdarnini 1984 




Davis et al. 1992b 




Taylor fc 


Rowan- 


Robinson 1992 




Klypin et al. 1994 




ling et al. 1994 


). 



At present this is a perfectly sa tisfactory solution (especially 
if a mo dest tilt is also allowed ( Liddle fc Lyth 1993b ; ISchae- 
fer fc ^afi 1994 jPogosyan fc Starobinsky 1995| ; [Liddle et 
al. 199q^ but loses the simplicity of having only cold dark 
matter. Even more fanciful solutions ex i st, including decay- 
ing neutrinos ( Bond fc Efstathiou 1991; Dodelson, Gyuk fc 



Turner 
cock ic|95| 
et al. 138 
Polarsk i 



1994 
)87|; 



White, Gelmini fc Silk 1995a 



McNally fc Pea 



broken power laws or double inflation (Turnei 



Gottlober, Miicket fc Starobinsky 1994 [Peter 



additional isocu rvature component (Kawasaki, Sugiyama fc 
Yanagida 1996; stompor, Banday fc Gorski 1996). There is 



Starobinsky 1994[ ; Kates et al. 199£ ), and an 



also the possibility that astrophysical processes could play 
a role i n determining the power spectrum of galaxy fluctu 



ations (Babul fc White 1991 



Efstathiou 1992 ; Bower et al 



1992 ; L^mbas et al. 1993[ ), although it is not at all clear that 



such effects will be significant 

A more fundamental change is to stay with CDM but 
abandon critical-density models, going instead to low den- 
sity. Viabl e models of this type are possible both i n open 
univers es (Ratra fc Peebles 1994 ; 36rski et al. 1995 ; Liddle 
et al. 1996a; Vamamoto fc Bunn 199t ) and in flat universes 



with a cosmological constant (Peebles 1984; Turner et al 



linear power spectrum (Klypin, Primack fc Holtzman 1996), 
continue to do so. 

So, although there are many alternative paths to ex- 
plore, it seems to us premature to relinquish the aesthetic 
appeal of critical density CDM models. Within the context 
of this general paradigm, there remain a number of param- 
eters to be honed by the data. The value of the baryon frac- 
tion was traditionally a flxed quantity, and also thought to 
be low enough that the cosmological effects would be rela- 
tively unimportant. This is no longer true if one allows for 
the higher values of SIb which now seem more credible, and 
moreover such variations are now significant at the level of 
detail which current data can distinguish. Raising Ob to 10 
per cent or above makes a substantial difference, as we will 
now show. 



3 THE DARK MATTER POWER SPECTRUM 

The matter power spectrum in the CDM model is spec- 
ified by two functions. The initial spectrum of perturba- 
tions is taken to be a power-law Pinit(fc) oc fc", where n 
is the slope. The case n = 1 is the scale-invariant spec- 
trum, and n 1 models are often referred to as 'tilted'. 
Infiationary models typically give n around one, though the 
possible range covers all the region of interest for struc- 
ture formation^. The transfer function T(fc) measures the 



t Inflationary models may also produce a spectrum of gravita- 
tional waves which can influence the COBE normalization; typ- 
ically these make it more difficult to fit the data and we do not 
consider them here. 
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Figure 1. Linear theory constraint curves for CDM models in the h~n plane for different values of Q,^ln? , namely (a) 0.0125, (b) 
0.025, (c) 0.0375 and (d) 0.05. The squares indicate where the models in Table 1 lie. The curves and shaded areas show 95 per cent 
confidence regions from the galaxy correlation function (solid), the height of the CMB peak (dashed) and the abundance of damped 
Lyman alpha systems (dotted). It is the tension between these constraints which is eased by allowing a higher fie- Cluster abundance, 
quasar abundance and velocity information are less constraining than these. The shaded region is allowed by all data constraints, and 
some continuation above the dashed line would be permitted if reionization was sufficiently early. 



amount of growth (relative to the infinite wavelength mode, 
and computed in linear theory) experienced by the modes of 
wavenumber k. The present linear power spectrum is given 
by P(k) oc Pinit(fc) r2(fc). In n = 1 CDM models, the trans- 
fer function is time-independent at late times; however, its 
form depends on the cosmological parameters h and Ob 
describing the matter content of the Universe. Commonly, 
CDM models are described via a fit to the standard forirj^ 
of Bardeen et al. (|l98e|) 



ln(l + 2.34g) 
2Mq 



(2) 



[l + 3.89g + (le.lg)'' + (5.46g)^ + (6.71g)^] 



-1/4 



with q = k/hV ( pugiyama 1995| ; |Hu fc Sugiyama 1996| ). The 
parameter F, a function of both h and JIb and referred to 
as the 'shape parameter', describes the shape of the matter 
transfer function. We have calculated transfer functions for 
our models numerically and then fit to this form. 

In observational papers, the best-fitting model for a par- 
ticular data set is often given in terms of an effective shape 
parameter. Toff. This is usually done by fitting to the power 

■t This form is not a good approximation at high unless 
h ~ 0.50. This could be important when considering ow Sin mod- 



els shou ld Vlvt turn out to be 10 per cent or higher (Blumonthal 
Dekel fc| Primack 1988|). 



spectrum under the assumption that n = 1. Since we will be 
allowing n different from unity, and also since the fit to the 
shape will depend on the range of scales considered, there is 
some ambiguity in the meaning of this Fcff , and so we will be 
explicit in our definition here. The 'effective' value of F for 
the theoretical models which we show in the table is obtained 
by finding the scale-invariant F-model which gives the same 
ratio of power on scales 8/i~^Mpc and 50/i~^Mpc, those 
scales being chosen as roughly the range typically probed 
by galaxy surveys. Since fitting a scale-invariant F-model to 
a tilted model is only approximate, it must be treated as 
only the roughest of guides and we shall always work from 
the true power spectrum. 



Following Liddle & Lyth (1993a), the normalization can 
be specified by considering the value of the density fluc- 
tuation at horizon-crossing: 5^{k) oc This is defined 
through the contribution to the rms fluctuation per loga- 
rithmic interval in wavenumber k as 



k^Pik) 
27r2 



(fe) T\k) , 



(3) 



which is dimensionless. The normalization can be speci- 
fied by giving the value of (5h at the present Hubble scale. 
The four- year COBE normalization of tilted critical-density 
CDM models can be well represented by a fitting function 



&(fc = i?o) 



(1.94 ±0.13) X 10" 



(4) 



© 1996 RAS, MNRAS 000, 000-000 



CDM models with high baryon content 5 



C 0.9 - 




0.05 



Figure 2. As Fig. 1, but in the Q,^-n plane with h fixed at 0.50. 
The small region to the bottom right is ruled out because it has 
insufficient power to produce damped Lyman alpha systems at 
2 = 4. The continuation of the shaded region upwards across the 
dashed line would be permitted by early reionization. 



exp [-0.95(n 



1) 



0.169(n 



which is accurate to b etter than 3 per cent for 0.7 < n < 1.2 
( |Bunn fc White 1996| ). The effects of VLb and h on this fit 
are negligible. 

A quantity related to the power spectrum is the variance 
of the density field smoothed on a scale 7?, which is simply 
the integral of A^, 



(5) 



where we take the smoothing function W{kR) to be a spher- 
ical top-hat, given by 



(6) 



Ways of comparing models with large-scale structure 
data are well developed in the literature. We use obser- 
vations as described in Liddle et al. (1996b,c) to outline 
the preferred regions of parameter space (see also Silk 1994; 
Bond 1994; Bahcall 1995; White et al. 1995b; Primack 1996). 
Four crucial types of data turn out to be the most con- 
straining, namely: the COBE observations of large-angle 
anisotropies; a compilation of intermediate-scale microwave 
anisotropy data; the abundance of damped Lyman alpha 
systems seen in quasar absorption spectra; and the shape 
of the galaxy correlation function. We shall consider several 
other types of data in addition to these, which prove not to 
give firm additional constraints. 

We shall describe all of these in more detail in the fol- 
lowing sections. However, to guide our later discussion, we 
show the results of the linear theory comparison now. The 
three parameters we can vary are Qb, h and n. Fig. |l| shows 
cuts of the parameter space at different choices of constant 
Qsh^- In each case there is an allowed region, but for the 
old nucleosynthesis value of f^B, shown in Fig. ^(a), agree- 
ment is available only for very low h. There is considerable 
pressure from direct measurement against reducing h below 
0.5, and it is rapidly problematic for large-scale structure 
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Table 1. Power spectrum measures for selected CDM models. 
The models have f2o = 1, h = 0.5, n = 0.8, with four values 
of fiB as shown. The values of F come from fitting to numerical 
evaluations of the transfer function for 10"'^ 1^ k < l/iMpc~^. 
The values of Fefi are fits to the ratio of the dispersions an at 
R = 50h~^ and R = 8h~^Mpc. The errors on aji reflect the 
1(T uncertainty from the COBE four-year normalization. The es- 
timates of the reionization redshift 2^; and optical depth t are 
highly uncertain. 



to go much above it, whatever the baryon density. For this 
reason we also show, in Fig. ^ a parameter-space cut at the 
value h = 0.5. This highlights the importance of varying the 
baryon density, and in many ways this figure can be said 
to represent our main result. Only by raising the baryon 
density to 10 per cent and above can critical-density CDM 
models be viable at /i = 0.5. Both the galaxy correlation 
function and the CMB peak height exert pressure in that 
direction. 

In order to focus our discussion, we select some partic- 
ular models for closer study. We choose to concentrate on 
h = 0.5, for which n — 0.8 is around the optimal choice, 
although less tilt can be accommodated for lower values of 
h. Our main aim is to consider variations of fls, and so we 
pick four values, corresponding to: the old BBN number (5 
per cent); a number consistent with the low deuterium mea- 
surements and more recent BBN calculations (10 per cent); 
the favoured number from the cluster baryon fraction (15 
per cent), which is some way above even the highest stan- 
dard nucleosynthesis upper limits; and an extreme example 
representing close to the top of the range of observed cluster 
baryon fractions (20 per cent). In Table |l| we show a sum- 
mary of important quantities for models spanning this range, 
calculated primarily in linear theory by numerical evolution 
of the coupled fiuid, Einstein and Boltzmann equations^ 
Here the quantity z^i is the redshift of reionization; it is es- 
timated using the techniques of Liddle & Lyth (1995), and is 
highly uncertain. The optical depth to Thomson scattering, 

r = 0.035 ^B ft [{l + z,if^^ 



1] 



(7) 



inherits this uncertainty, and measures the probability of mi- 
crowave photons being re-scattered, assuming full ionization 
from 2ri to the present. 



3 In White et al. (1995b) the values for ai were obtained from 
extrapolating P{k) as a power law from 8/i~'^Mpc down to 
lh~^Mpc. This neglects the curvature of the transfer function 
over this range of scales and over-estimates tri. We use the full 
transfer function here. 
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Figure 3. The matter power spectrum for £7b = 0.05,0.1,0.2, 
the highest lines having the lowest baryon density (we omit the 
15 per cent model for the sake of clarity). All models have /i = 0.5 
and n = 0.8. The solid squares are taken from Peacock & Dodds 
(1994) with the 4 points at highest k omitted. We have adjusted 
the data to the best fit normalization for the =0.1 model. 
The open circle is the velocity-derived point from Kolatt & Dekel 
(1996); note however that the error bar shown on this point omits 
cosmic variance, which is large. 



Figure 4. The redshift space matter power spectrum for the 
fie =0.1 model from Table 1 (solid line). The dashed line shows 
the result for standard CDM. The data are taken, without ad- 
justment, from Tadros & Efstathiou (1995; weight 8000); approx- 
imately every second point is uncorrelated. This ignores a down- 
ward correction to the theoretical curves of roughly 15 per cent 
for k < O.l/iMpc"^, arising from the window function of the 
surveys. 



4 STRUCTURE FORMATION 

We now discuss the main limits imposed on cosmological 
models from observations of the clustering properties of 
galaxies, the galaxy cluster abundance, peculiar velocities 
and high redshift object formation. 



4.1 The galaxy correlation function 

Perhaps the classic problem of standard CDM has been 
that the model has the wrong 'shape'. Specifically the ra- 
tio of large-scale (~ 50/i~^Mpc) power to small-scale (~ 
8/i"^Mpc) power is too small compared to observations. 
Both tilt and small-scale damping from a high baryon frac- 
tion help address this. 

The current situation regarding the data relating to 
shape is somewhat uncertain. If we co mpar e our models to 
the compilation of Peacock & Dodds (1994), taking the er- 



ror bar* hi fade value, thfetl wte I'tequite S7b > 0.1 tO prftvidt; 
sufBcieiiL laige-scale puwei, as can be seen in Fig. 3. Fbi 
the model with ~ 0. 1 and n = 0.8, a simple fit to the 



Peacock & Dodds (1994) data points assuming Gaussian un- 
correlated error bars and allowing the normalization to vary, 
shows the model is excluded at about 97 per cent confidence, 
if one drops the four shortes t-scale points which may be af- 
fected by non- linear biasing ( Peacock 199(: ). If all the points 
are included this model fits at the 95 per cent confidence 
level. The models with higher fie fare better, being allowed 
at better than the 95 per cent confidence level, even exclud- 
ing the four shortest-scale points. A reasonable conclusion 
based on this constraint then is that critical-density CDM 
models (with h ~ 0.5 and n ~ 0.8) are allowed provided 
that JIb is about 10 per cent or gr eater. 

The Peacock & Dodds (1994) data are a compilation 
to which several adjustments and/or corrections have been 



made. Since we have specific models in mind, we can carry 
the predictions forward to confront individual data-sets di- 
rectly, rather than comparing to processed data. An example 
is shown in Fig. ^, where we compare the Ob =0.1 model 
to the power spectrum of IRAS galaxies in redshift space, 
conveniently tabulated in Tadros & Efstathiou (1995). We 
have chosen this fls as the model which is only marginally 
allowed by Peacock & Dodds ( |l994| ) . In this model the IRAS 
galaxies are almost unbiased tracers of the mass - we have 
not scaled the theory or the data points in Fig. ^ in any 
way. Specifically, we have assumed a redshift-space distor- 
tion parameter /3 = 1, and a small-scale exponential velocity 
distribution with a = 280 km s^^ to convert from real to red- 



shift space (Cole et al. 1995). It is clear from Fig. ^ that this 
model fares extremely well against this survey, though it is 
less constraining than that of Peacock & Dodds (1994) since 
it incorporates less data. 

The first firm indication that standard CDM lacked 



power on large scales came from the APM survey ( Mad 



dox et al. 1990| ; [Maddox, Efstathiou fc Sutherland 1996| )."Tb 
compare to this survey, we show in Fig. ^ the predicted an- 
gular correlation function w( 8) for our four models. We use 
the fit in Jain, Mo & White (1995) to obtain the non-linear 



power spectrum; this is particularly important on scales less 
than 1°, but also has a modest effect at larger angles where 
it suppresses the correlation slightly. All our models clearly 
have sufficient large-scale power to fit this data, in contrast 
to the standard CDM model. 



4.2 Galaxy cluster abundance 

The present abundance of galaxy clusters requires ag — 0.6, 
with most authors agreeing fairly well on the centr al value 
but with differing opinions as to the uncertainty (Evrard 



198S|; [Henry fc Arnaud 199l[ [Bahcall fc Cen 1993[ [White 
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Figure 5. The galaxy angular correlation function w{9) for the 
APM survey as predicted by our models with h = 0.5, n = 0.8 
and Hb = 0.05 (long dashed), 0.10 (solid), 0.15 (dotted) and 0.20 
(short dashed) . All curves have been COBE normalized with bias 
chosen to agree with the data at ~ 1°. They include a correc- 
tion from the non-linear part of the power spectrum. The shaded 
region represents the 2a error band from a fit to the APM data 
(Baugh & Efstathiou 1994). 



et al. 1993a; 



Bond fc Myers 1996c; Viana fc Liddle 1996 



Eke, Cole & Frenk 1996). Many determinations do not per- 
mit values as high as our lowest-b aryon model predicts (see 
Table |l|), though Viana & Liddle ( |l996| ) quote -f 32 per cent 
and -24 per cent at 95 per cent confidence, among the most 
conservative, which does just allow it. However, things are 
much better for the models with a higher baryon content. 
We find that any model of the type considered in this paper 
satisfying both the COBE and the galaxy correlation data 
constraints automatically satisfies the cluster constraint as 
well. 



4.3 Peculiar velocities 

The galaxy peculiar veloci ty field can be recovered th rough 
the POTENT procedure ( [Bertschinger fc Dekel 198£| ) from 
pecuUar velocity catalogs, such as the Mark III catalogue. 
It can then be used to constrain the amplitude of the power 
spectrum either directly, th rough the spatial derivatives of 
the p eculiar velocity field (Dekel 1994; Strauss fc Willick 
1995|), or iri directlv, by means of the bulk velocity " ^Kok^T^ 



Dekel 1996 ). Because the latter is sensitive to larger scales. 



results obtained through both methods do not necessarily 
coincide. Due to the very large cosmic variance, resulting 
from the fact that only local measurements have been made 
from a random field, we find that neither method further 
constrains the available parameter space left by the other 
observations. However, for comparison, we show in Fig. ^ a 
point fro m the POTENT velocity power spectrum of Kolatt 
& Dekel (1996). Notice that our models fit this well even 
without adding the cosmic variance. We note in passing that 
the velocity data can also be used to constrain the shape of 
the power spectrum, though at the moment results are not 
as robust as those using galaxy correlation data. Zaroubi et 
al. (1996) seem to prefer a higher Fcff , between 0.35 and 0.55 



when marginalized over era, than the galaxy surveys, though 
the uncertainty is large. 

In addition, the question of the small-scale ('^ 
lft~^Mpc) velocities should be addressed, i.e. is the veloc- 
ity field cold enough? This problem has been discussed by 
many people (see White et al. 1995b for references) but a 
clean compar ison w ith the data remains difficult. Following 
White et al. (1995b), we compute the second moment of the 
mass correlation function J2(l /i~^Mpc), which is related to 
the rms velocity on the same scale. We find a reduction in 
short-scale power, as measured by J2(l /i^^Mpc), by a factor 
of 4 to 5 for our JIb 0.1 models, compared with standard 
CDM. Although this is a significant reduction, it may still 
be too hot relative to that measured in redshift catalogues, 
though this may be subject to substantial cosmic variance. 
We believe that the situation regarding velocities on these 
scales will probably require careful comparison of A''-body 
simulations with observational results, which is beyond the 
scope of this paper. 



4.4 High-redshift object abundance 

We are going to concentrate on three types of objects found 
at high redshifts: damped Lyman alpha system (DLAS), 
galaxies and quasars. We will also mention the situation 
regarding clusters at intermediate redshifts. The technique 
we use to compare observed abundances with the one s pre - 
dicted by our models is the standard Press-Schechter (1974) 
calculation (see Liddle at al. 1996b for a discussion of our 
particular implementation) . 



4-4- i Damped Lyman alpha systems 



The observations of DLAS by Storrie-Lombardi et al. (1991: 



imply that at redshift 4, the most constraining point, the 
amount of neutral hydrogen in bound objects is 



fim = (0.0011 ±0.0002) h~ 



(8) 



To obtain a conservative bound, we assume that all the hy- 
drogen in DLAS is in the neutral state, and taking the 2a 
lower limit, we obtain a lower bound on the baryon density 
at that redshift, 



^Biz = 4) > 0.0007/1" 



(9) 



In order to compare this limit with our models we need to 
make some further assumptions regarding the nature of the 
DLAS. To be conservative we assume that the DLAS at 
such a high redshift arise predominantly within filamentary 
structures, which see ms to be support ed by recent hydrody- 
namical simulations ( Katz et al. 1996| ) , and consequently we 
use the threshold value of Sc — 1.5 in the Press-Schechter 



calculation ( Monaco 1995 ). We also conservatively use the 
95 per cent upper limit on the COBE normalization, though 
the models also fit well at the central normalization. 

In Fig. ^ we show an estimate of the baryon density in 
objects with a mass in excess of M at redshift 4 for our mod- 
els. There is little difference at high redshifts between the 
models with higher baryon fractions. These models would 
come under pressure were it shown that the observed sys- 
tems correspond to substantially more massive objects than 
usually supposed. However, under the reasonable assump- 
tion that the DLAS correspond to objects of masses 10^" 
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Figure 6. The baryon density in objects of mass greater than 
M at redshift 4, for our models. Along the right-hand side, the 
higher lines correspond to lower baryon density, the lines being 5 
(long dashed), 10 (solid), 15 (dotted), and 20 (short dashed) per 
cent. The horizontal dotted line indicates the 95 per cent lower 
limit from DLAS observations (Storrie-Lombardi et al. 1995). 



Mq and greater, all of the models are comfortably al- 
lowed. We show the full parameter space constraints as the 
dotted curves in Figs. |l] and The DLAS provide a limit 
on how much tilt can be introduced; at high baryon density 
this is presently a stronger constraint than the CMB peak 
height, as seen in Fig. |^. 



Lyman break galaxies 

Recently Steidel and collaborators (Steidel at al. 1996a,b; 
Giavalisco, Steidel & Macchetto 1996) were able to detect 
galaxies between redshifts 3 and 3.5 using a technique based 
on the identification of the redshifted Lyman continuum 
break. Their results put a lower bound on the number den- 
sity of Lyman break galaxies in that redshift range: 

iVgai(2 = 3.25) > 0.00288 Mpc"^ . (10) 

In Fig. 1^ we plot the number density of collapsed objects 
with mass above M at redshift 3.25, for our four specific 
models. Here we use the threshold value of 5c = 1.6 in the 
Press-Schechter calculation, which seems the optimal when 
the objects under consideration might have collapsed pre- 
dominantly non-spherically due to local tidal fields (Monaco 
1995; Bond & Myers 1996a,b). As before, we conservatively 
use the 95 per cent upper limit on the COBE normalization. 

In order to constrain our models we need a rough es- 
timate of the mass of the virialised halos with which the 
Lyman break galaxies are associated. According to the ob- 
servations (Steidel et al. 1996b, Giavalisco et al. 1996) the lu- 
minosity of these galaxies arises predominantly within an ap- 
proximately spherical central region with a radius of about 
2.5 h^^ kpc, where the velocity dispersion is between 180 to 
320 kms"'^. Assuming that galaxy halos can be roughly ap- 
proximated by a truncated singular isothermal sphere, these 
velocity dispersions correspond to halo masses in the region 
of 5 X 10" /i"^ Mq to 3 X 10^^ Mq. However, the quoted 
interval for the velocity dispersion was obtained assuming 



Figure 7. The number density of collapsed objects with mass 
above M for the four models at redshift 3.25. Higher lines corre- 
spond to lower baryon density, the lines being 5, 10, 15 and 20 
per cent. The dotted line indicates the observed number density 
of Lyman break galaxies (Steidel et al. 1996a). 



that the observed broadening of interstellar absorption lines 
is solely due to gravitational induced motions. If part of 
the broadening is due to non-gravitational effects, such as 
interstellar shocks arising from local star-forming regions, 
then the velocity dispersions would be smaller and hence 
also the derived halo masses. On the other hand, the true 
mean velocity dispersion of the halo is probably higher than 
the velocity dispersion of the stars, as these feel the effect of 
only the central part of the halo mass. In the future, spec- 
tra with higher signal to noise ratios and better resolution 
may enable the distinction between different sources of line 
broadening. 

Another pro blem fo r our models may be the observation 
by Steidel et al. (1996a), using additional K-band photome- 
try, that the Lyman break galaxies appear redder than one 
would expect if they were undergoing their first burst of 
star formation. However, there remains considerable uncer- 
tainty in estimating formation time from colours. Neverthe- 
less, were it shown that these galaxies formed substantially 
earlier than the redshift at which they are observed, all mod- 
els with critical density would be under severe strain. 



4.4-3 Quasars 

For a long time, the abundance of quasars was the only in- 
formation available about the process of structure formation 
at high redshifts. Now the abundances of DLAS, and more 
recently of Lyman break galaxies, probe similar epochs. Fur- 
ther, it has been shown that in general both DLAS abun- 
dan ces (see e.g. Liddle et a l. 1996b) and Lyman break galax- 
ies (Mo & Fukugita 1996) provide stronger constraints on 
theoretical models of structure formation than quasar abun- 
dances. The observed quasar number density at a redshift 
of 3.25 is about 10"^ Mpc"^ (Schmidt, Schneider & 
Gunn 1995), and as can be seen in Fig. |^ this constraint 
is much weaker than that arising from the abundance of Ly- 
man break galaxies; all models are viable for any quasar host 
galaxy mass up to 10^'^ ^0- even higher redshifts 
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the observed quasar number density does become more con- 
straining, but not by much. For example, at redshift 4 our 
models are still able to reproduce the observed quasar num - 
ber density, about 3 x 10~^/i^Mpc"^ ( Schmidt et al. 1995 ), 
for any host galaxy mass up to 10^^ M0, except the 



model vith Q,b = 0.2 which requires host galaxy masses 
from 6 X 10^^ Mq upwards. 

We also carried o ut a m ore careful analysis along the 
lines of Nusser & Silk (1993), allowing for relations between 



the abundance of quasars and their mean lifetime, and the 
luminosity and mean lifetime of a quasar and its minimum 
mass. We found that the conclusions given above hold for 
quasar mean lifetimes over a broad interval between about 
lO'^ and 10^ yrs, the precise range depending on the fraction 
of collapsed halos hosting quasars, the amount of baryons in 
the halo used to form the central black hole and the quasar 
radiative efficiency. Our models and conclusions are similar 
to the tilted cases considered by Haehnelt (1993) and by 
Nusser & Silk (1993), who found that n ~ 0.8 models were 
allowed, although the quasars were then required to be ra- 
diating near to the Eddington limit and with lifetimes of 
~ 10^ years. Note that for quasars, in common with other 
high-z objects, there is a great degree of uncertainty in the 
estimation of the amount of collapsed matter at early times 
(see e.g. Katz et al. 1994; Eisenstein & Loeb 1995), and 
that it may be easier to form central black holes in higher 
JIb models. 



4-4-4 Clusters 

Finally, regarding intermediate redshift X-ray clusters, the 
best available observations at this moment are those of Lup- 
pino & Gioia (1995). Their sample of 6 clusters was selected 
from the Einstein Observatory Extended Medium Sensitiv- 
ity Survey (EMSS). They have luminosities in excess of 
Lx = 1.25 X 10'^'^h~^ ergs s~^ and redshifts between 0.55 
and 0.83, where the most distant ever detected X-ray se- 
lected cluster is located. Given the comoving volume over 
which the survey was carried out, they conclude that the 
comoving number density of X-ray clusters with luminosi- 



ties in excess of Lx is (8.8±3.44) x 10" ft Mpc""^, at a mean 
redshift of 0.66. 

Again the difficulty is in obtaining a reliable mass es- 
timate for these objects. However if we instead concentrate 
on the cluster X-ray temperatures, we can get a rough idea 
of what to expect from our models at 2 = 0.66. For that 
we first need to relate the X-ray luminosity Lx to an X-ray 
temperature. To a first approximation we can assume that 
the L — T relation at z ero re dshift also applies at z = 0.66. 
Following Henry et al. (1992), we translate Lx into a bolo- 



metric luminosity using Lboi ~ 1.18(fcr)'''^^I/x, where kT is 
in keV. Using the observed zero redshift L~T relation given 
in David et al. (1993), and assuming h ~ 0.5, we find that Lx 
corresponds to an X-ra y tem perature of roughly 6keV. Fol- 
lowing Viana & Liddle (199(: ), and using the values of as in 
Table ^, we calculate the mass corresponding to these clus- 
ters under the constraint that their observed zero re dshift 
number density be recovered (Henry & Arnaud 1991). We 



assume Sc ~ 1.6 in the Press-Schechter calculation, integrate 
over temperature and obtain an estimate of the number den- 
sity of X-ray clusters with Lboi > 1.1 x 10*^/i~^ergs s~^ at 
z — 0.66. For our models with Ob ~ 0.05 and f^B ~ 0.20, we 



obtain 4 x 10*^ and 2 x 10*^ x 10"^/i^Mpc"^ respectively, 
with the other models lying in between. Here the uncer- 
tainty, an order of magnitude in each direction, comes from 
pushing both the COBE normalization and the observed 
present numl jer density of 6keV X-ray clusters ( Henry & 



Arnaud 1991) to their 2a confidence limits. As can be seen. 



all our models have mild difficulty in producing the observed 
number density of high-redshift clusters, but in view of the 
uncertainties there is still plenty of room for manoeuvre. 
However, current data may suggest that at fixed tempera- 
ture the X-ray cluster luminosity is smaller at higher red- 
shifts (e.g. Navarro, Frenk & White 1995) ; if this trend turns 
out to be correct, then all models with critical density will 
have great difficulty in reproducing the observed number 
density of high-redshift clusters. 



5 COSMIC MICROWAVE BACKGROUND 
ANISOTROPIES 

Any theory which purports to explain structure formation 
in the universe must fit both the large-scale structure data 
and the data on CMB anisotropies. The most influential 
piece of CMB data is the COBE normalization of large-scale 
structure theories (see White & Scott (1996) for a discussion 
and list of references); however, with the influx of data on 
smaller angular scales there are now extra constraints which 
must be met. For our purposes the strongest of these is that 
our model must be able to produce sufficient degree-scale 
power, which limits the amount to which we can tilt the 
primordial spectrum to alleviate the large-scale structure 
problems. 

We show the radiation anisotropy power spectrum for 
our models, normalized to the COBE four-year data, in 
Fig. 1^, along with a standard CDM model for comparison. 
The angular power spectrum is plotted in the usual way, 
with £{£ + l)Ce in dimensionless units. The index of the 
multipole expansion £ ~ , with 6 in radians. We have 
chosen not to include any late reionization in these mod- 
els, as the estimate of the reionization redshift is very un- 
certain. Including reionization would reduce the power on 
scales £ 100 by exp(— 2r), where estimates of r are quoted 
in Table 1. 

Notice that the increase in the baryon content has mod- 
ulated the height of the peaks in the spectrum. The flrst and 
third peaks are increased in amplitude while the second peak 
is decreased. The tilt of the spectrum has monotonically de- 
creased the power at small angular scales (high £). This lack 
of power on smaller angular scales, ^ 15', is a robust fea- 
ture of our models which should be testable by the next 
generation of interferometers and array receivers. 

Also shown in Fig. ^ is a selection of the current data 
on CMB anisotropies. From left to right the solid circles are: 
COBE (Hinshaw et al 1996; 3 points, 1 upper lim it), SP94 
( iGundersen et al. 1995| ), MAX ( Tanakaet al. 1996| ), Python 
(Ruhl et al. 1996, 3-beam), MS AM (Inman et al. 1996, 3- 
beam) and CAT (Scott et al. 1996, 2 points). The error bars 
arc ±1(7 while the horizontal lines indicate the range of an- 
gular scales to which the experiment is pr imarily sensitive 



The solid squares indicate the SK95 data (Netterfield et al 



1996), with no allowance made for the 14 per cent calibration 
uncertainty. The upper limits are the highest-^ point from 
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Figure 8. The radiation power spectrum for our models (dashed) with f2o = 1, h = 0.5 and n = 0.8. The curves have Qb = 0.05, 0.1 
(thick) and 0.2 from bottom to top at £ ~ 200 (we have omitted the Hb = 0.15 Une for clarity). The solid line is sCDM, with n = 1, 
h = 0.5 and 5 per cent baryons. All models have been normalized to the 4-year COBE data. The points are recent la detections and 
95 per cent upper limits as discussed in the text. Solid circles from left to right are: COBE (3 points, 1 upper limit), SP94, MAX, 
Python, MSAM, CAT (2 points). White Dish, OVRO, ATCA. The solid squares are the SK95 points, without inclusion of the calibration 
uncertainty in the error bars. The error bars shown on this figure should be regarded as indicative only. 



COBE at ^ ~ 30, t hen at large £ are Wh ite Dish ( [Tucker ct 
al. 199 J), OVRO (iReadhead et al. 1989|) and ATCA (^ub 



rahmar yan et al. 1993 ), all quoted at 95 per cent confidence. 
For display purposes only, the points have been converted 
from temperature to power (temperature squared) with er- 
ror bars (Ty2 = 2T(jt- 

These data have been selected from observations which 
have shown repeatability and which probe angular scales 
near the first peak. We have not included any allowance 
in the error bars for the effects of foreground subtraction. 
The expected or measured size of the foregrounds is dis- 
cussed in the papers quoted above. Typically the experi- 
ments ba-^re mapped regions for which foreground contami — 



nation k known tn hp low and /or have used tipvpral frpquen- 

cies to test for the presence of foregrounds. The subtrac- 
tion or correction for foregrounds will generally weaken the 
lower limit on n. Where experiments have quoted two num- 
bers from the same patch of sky we have chosen the high- 
est I, which provides the largest lever arm. This is because 
such points are correlated in an unknown way. However, due 
to the observing geometry and data reduction st rategy, the 
SK95 p oint s are less than 20 per cent correlated (Netterfield 
et al. 1996), except for a common 14 per cent calibration 



uncertainty, so we have included all 5 as uncorrelated. Cor- 
rections to this assumption should also permit a lower n. 

To test the constraining power of the CMB data on our 
theory, we have computed the likelihood of fitting the data 
as a function of overall normalization (Q), tilt n and an addi- 
tional multiplicative factor for the SK95 points to represent 
the calibration uncertainty. We have not tried to include 
the upper limits or the CAT points (which are correlated 
and have very large error bars). We show the likelihood in 
Fig. ^, marginalized over the normalization and SK95 cal- 
ibration uncertainty. We have chosen a flat prior for (Q), 
which is anyway well determined by the COBE data. For 
the calibration uncertainty of SK95 we multiply by a Gaus- 
sian centre d on unity with width 14 per cent (Netterfield 



et al. 1996). The calibration uncertainty from the other ex- 



periments was added in quadrature to the error bars on the 
points. 

Our analysis here has not been very sophisticated, but 
it shows that while the tilted model with 5 per cent baryons 
does indeed have d ifficulty fitting the height of the peak 
(White et al. 1995b), models with 10 to 20 per cent baryons 
fare quite well. It would be interesting to perform a fit of 
these models to all of the data sets, rather than just the 
band-powers, to see how well they fare in detail. The best 
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Figure 9. The marginal likelihood versus n, integrated over nor- 
malization and SK95 calibration, for a fit to a selection of the 
current CMB data. We show likelihoods (from right to left) for 
Ub = 0.05, 0.1 (solid), 0.15 and 0.2. In all cases, h is fixed at 0.50. 



fitting n depends on the assumptions made about the tails of 
the SK95 calibration uncertainty since tlie fit of any model 
to the data becomes much better as the SK95 data are 'low- 
ered'. The of the best-fitting model decreases by ~ 10 if 
the SK95 calibration is decreased from 1 to 0.82 (consistent 
with the MSAM calibration) for example. With the SK95 
cahbration fixed at 1 the best-fitting model fits at about 95 
per cent confidence level, this becomes 50 per cent if the 
calibration is fixed at 0.82. The first and second moments of 
the marginal likelihoods from Fig. |^ are n = 0.98 ± 0.06, 
0.92 ± 0.06, 0.86 ± 0.06 and 0.80 ± 0.06 for Qb = 0.05, 
0.10, 0.15 and 0.20 respectively. If one includes the possi- 
bility of reionization, then the n = 0.8 model with 10 per 
cent baryons begins to become unlikely. Including optical 
depth r, as given by Table ^, is approximately equivalent 
to shifting n — > n — 0.65r. Should observations further in- 
crease the preferred peak height, then the model with 10 
per cent baryons will be disfavoured. On the other hand, if 
the height were to decrease, e.g. if some of the anisotropy 
were non-cosmological, even more tilt or reionization could 
be tolerated. 



6 CONCLUSIONS 

In this paper we have investigated the effects of assuming 
a baryon fraction two or more times higher than 1993 esti- 
mates of BBN would predict. With Qb ^ 0.1 the baryons 
are becoming dynamically relevant, and so the increase has 
a non-negligible impact on the model of structure formation. 
We have concentrated on cold dark matter models with crit- 
ical density, allowing Qsh^ to vary within the range 0.0125 
to 0.05. 

Inclusion of a higher baryon fraction damps small-scale 
power in these models, since the growth in the baryons is 



coming into conflict with the growing evidence for enhanced 
power in the CMB on degree scales. 

In order to facilitate more complete comparison of our 
model with data and to encourage numerical experiments, 
we have chosen a 'preferred' model with h = 0.5, fls = 0.12 
and n = 0.8. This model is not the absolute best-fit when 
considering large-scale structure data alone, which would 
favour some combination of a slightly lower h, lower n or 
higher Q,b, as seen in Figs. |^ and|^. In selecting it, we have 
also paid heed to other evidence, such as nucleosynthesis and 
direct measurement of the Hubble constant, which suggest 
that one should avoid straying too far from the 'standard' 
values. Clearly changes in these values of Hb, h and n can 
play off against each other. Higher Ob allows less tilt, and 
vice versa. If reionization occurs early, or if future CMB data 
require more power on degree scales, then this limits the tilt 
and thus the highest h which can be tolerated. 

Several developments could falsify this model. Firstly 
we require that the Hubble constant be close to 50 kms~^ 
Mpc~^. Should it be higher than 55 kms~^ Mpc~^, it would 
be very difficult to fit the current data while retaining a criti- 
cal matter density and pure CDM. A higher Hubble constant 
is in any case difficult for a critical-density universe since it 
predicts an age of below 12 Gyr. Of course, any firm mea- 
surement of the total density fio which obtained less than 
one would also rule this model out. In a similar vein, our 
model predicts ag to be in the interval from 0.5 to 0.8, and 
hence that IRAS galaxies are almost unbiased tracers of the 
mass. Both results ought to be refutable using future analy- 
ses of, respectively, velocity fiows and redshift space distor- 
tions in galaxy surveys. Furthermore, the standard theory of 
big bang nucleosynthesis is only consistent with our assumed 
baryon fraction if the helium mass fraction Yp has been un- 
derestimated. Should a firm upper limit on Yp around 24.5 
per cent be made then this model would require a violation 
of standard homogeneous nucleosynthesis by the introduc- 
tion of extra physics a.t t ^ 3 minutes. We also require the 
low deuterium determinations to be essentially correct, with 
the implication of little destruction of primordial deuterium. 

There are several specific predictions which the high f^B 
critical-density model makes. Starting with the theoretically 
cleanest of these, we predict that the radiation power spec- 
trum will have reduced power, compared to standard CDM, 
at £ ~ 600, corresponding to angular scales ~ 15'. This is 
seen clearly in Figure ^ and should be probed by the new 
generation of interferometer experiments CBI, VCA and 
VSA, and ultimately by the MAP and COBRAS/SAMBA 
satellites. On large scales, the shape of the matter power 
spectrum will differ from a scale-invariant F-model, due to 
the tilt in the initial spectrum of perturbations. Since our 
model has critical density, objects do not form as early as 
in open or cosmological constant dominated cosmologies, 
though the redshift of object formation is subject to large 
theoretical uncertainties. Finally, we predict a much higher 
fraction of dark baryons, though we are unable to say what 
form they will take. However, it is tempting to imagine that 
the halos of galaxies may have around a 20 per cent contri- 
bution from low-mass stars and other compact objects. 



the recent results of the MACHO experiment indicate ( Ben- 



suppressed by their interactions with the (relativistic) pho- 



nett et al. 1996b). We would certainly find it surprising if 



tons. More importantly, the higher baryon fraction allows 
one to tilt the initial spectrum by a greater amount before 



the MACHO fraction in galaxies was significantly less than 
the baryonic fraction in the Universe as a whole. 
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We briefly comment on the effect an increased baryon 
density would have on other structure- formation models. In 
all cases, the effect of the change to the CMB peak height is 
much more important than the change to the matter power 
spectrum. Recall that for critical-density CDM models, the 
main benefit is an increased CMB peak height, which is nec- 
essary to compensate for the need to tilt in order to get the 
galaxy correlations right. In low-density models, the spectral 
shape is already corrected by the changed value of floh, and 
so tilt is not required. However, observations only weakly 
constrain the tilt in such models, and so changing the baryon 
density will have a fairly neutral effect, offering neither ad- 
vantages nor disadvantages — it simply moves the allowed 
region of parameters around a little. For high fls/i^o, it may 
be possible to see oscillations in the matter power spectrum 
from the baryon component, but at critical density these are 
extremely small even in the 20 per cent model. For cold plus 
hot dark matter models on the other hand, a high baryon 
density could be threatening, because one must have less 
tilt in those models (the hot component having already ad- 
justed the spectral shape) and hence there is the danger of 
too high a peak in the CMB. Early reionization might then 
be necessary to reduce it to acceptable levels. 

Finally, let us end by noting that although critical- 
density CDM models remain viable, the present observa- 
tional situation allows only a very limited range of param- 
eters. If this paradigm is to be believed, then already the 
value of h is fixed at about the ten per cent level near 0.50, 
and similarly the slope n has to be around 0.8 ± 0.1. 
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